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1. INTRODUCTION 

In the standard hot big bang model, relics from the early universe can give us much 
useful information about microphysics. For instance the abundances of the light elements, 
produced when the universe was at a temperature of ~ 1 MeV, told us long before the 
existence of LEP that there were at most four, and probably three, species of neutrinos 
[Ij]. The most obvious of big bang relics are the baryons, which make our own existence 
possible. Furthermore, the universe seems to contain relatively few antibaryons. There 
is clear evidence that at least the local cluster of galaxies is made of matter, and there 
is no plausible mechanism to separate matter from antimatter on such large scales. The 
observed abundance of baryons today implies that when the universe was much hotter 
than a GeV the ratio of antibaryons to baryons must have been about one part in 10 8 

Sakharov pointed out in 1967 that this cosmological asymmetry between matter 
and antimatter could be a calculable result of particle interactions in the early universe, 
teaching us several profound things about fundamental physics: 

1) The number of baryons is not conserved. 

2) C (charge conjugation symmetry) and CP (the product of charge conjugation sym- 
metry and parity) are not exact symmetries. 

3) The universe must have been out of thermal equilibrium in order to produce net 
baryon number, since the numbers of baryons and of antibaryons are equal in thermal 
equilibrium. Note that all known interactions are in thermal equilibrium when the 
temperature of the universe is between 100 GeV and 10 12 GeV. 

Sakharov suggested that baryogenesis took place immediately after the big bang, at a tem- 
perature not far below the Planck scale of 10 19 GeV, when the universe was expanding so 
rapidly that many processes were out of thermal equilibrium. This scenario was explic- 
itly realized with the advent of grand unified theories (GUTs) 0-0, which predict baryon 
number violation, and possibly CP violation, mediated by ultra-heavy particles |L0[ . 

This "standard" GUT baryogenesis scenario is appealingly simple but does not easily 
fit into an acceptable cosmology. One difficulty is that in the standard model of electroweak 
interactions baryon number is known theoretically to be anomalous and not exactly con- 
served [|11[| . At low temperature this anomalous baryon violation only proceeds via an 
exponentially suppressed tunnelling process, at a rate proportional to exp(— An / a w ) ~ 0. 
At temperatures above ~ 100 GeV, however, electroweak baryon violation may proceed 
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rapidly enough to equilibrate to zero any baryons produced by grand unification |12|j . 
In the last decade there has been much work indicating that this is the case [|T3| - pO| . How- 
ever anomalous electroweak processes conserve B — L, the difference between baryon and 
lepton numbers, and so a net B — L generated by grand unified or other interactions will 
not be erased by electroweak interactions. In the minimal standard model if B — L is non 
zero then the equilibrium baryon number at high temperature is B = (12/37) (B — L) [fl5| , 
The possibility that the observed baryons are a relic from a net B — L generated 



at the GUT scale would be constrained if new B — L violating physics (such as Majorana 



neutrino masses) were discovered 24, E5l p7H29|. 



Another problem with GUTs is that they predict a relic abundance of massive stable 
magnetic monopoles much larger than the observed matter density [|3~0 , 31 1. The simplest 
GUT models which have an experimentally acceptable weak mixing angle are supersym- 



metric [32], |33[], and in these models the relic abundance of gravitinos (supersymmetric 
partners of the graviton) is incompatible with constraints on the mass density of the uni- 
verse f34T-|39|j . The most attractive solution to these problems, which also explains the 



flatness and homogeneity of the observed universe, is inflation [^0| - ^2| , which greatly di- 
lutes the abundances of unwanted relics. However after an inflationary epoch any remaining 
asymmetry between the numbers of baryon and antibaryons is negligible, and the universe 
typically reheats to a temperature which is less than or of order 10 12 GeV, well below the 
scale of grand unification |]J . Furthermore gravitino decays may affect nucleosynthesis [p7| , 
38fl unless the reheat temperature is below 10 9 GeV, a temperature too cool to reinstate 
the baryon asymmetry through GUT processes. 

There have also been many explanations of the matter-antimatter asymmetry which 
do not rely on GUTs |[26|| . Typically these propose the existence of new interactions and 
particles which exist solely in order to create baryons, and which are otherwise unmotivated 
and experimentally untestable. 

In 1985 Kuzmin, Rubakov & Shaposhnikov suggested an elegant and simple solution 
to the baryogenesis problem [T5|] . They argued that anomalous baryon number violation in 
the standard model of electroweak interactions is rapid at high temperatures and that the 



weak phase transition |43| , if first order with supercooling, provides a natural way for the 
universe to get out of thermal equilibrium at weak scale temperatures. Eventually bubbles 



1 There has been some confusion in the literature over this calculation; the correct fraction is 



found in [25, 26|. 
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of the broken phase nucleate and expand until they fill the universe; local departure from 
thermal equilibrium takes place in the vicinity of the expanding bubble walls. Since C and 
CP are known to be violated by the electroweak interactions, it is possible to satisfy all the 
Sakharov baryogenesis conditions within the standard model. This suggestion opened up 
the exciting possibility that we may be able to compute the baryon number of the universe 
in terms of experimentally measurable aspects of the weak interactions, such as particle 
masses and CP violating parameters. 

The effects of the CP violation in the minimal standard model from the Kobayashi- 



Maskawa (KM) phase |44| are much too small to explain the observed baryon to entropy 
ratioil. Explaining this number via baryogenesis in the early universe would require a C 
and CP violating asymmetry larger than 10 -8 in some out of equilibrium reaction, while 
at high temperature CP violation from the KM phase is always multiplied by a function of 
small couplings and mixing angles of order 10 -20 [45, 47]. For example, in minimal SU(5) 
GUT-scale baryogenesis with only KM CP violation one finds a baryon to photon ratio 
n B jn 1 ~ 10 -20 [[|. However the feasibility of baryogenesis during the electroweak phase 
transition has been demonstrated in several simple non-minimal versions of the standard 
model [p8|-p5|. Furthermore the baryon asymmetry of the early universe can be used 



to constrain parameters in such non-minimal models; for instance in the singlet majoron 



model ||56|| an acceptable baryon abundance is obtainable only if the mass of the r neutrino 
is greater than 5 MeV |^8[ . 

A crucial test of weak scale baryogenesis comes from the requirement that any baryons 
produced during the transition should survive until the present. If anomalous baryon vi- 
olating processes occur at a rate which is faster than the expansion rate of the universe, 
then thermal and chemical equilibrium will be restored after the transition. Since the 
equilibrium abundances of baryons and antibaryons are equal, any baryon number which 
was created during the transition could be washed out. In the broken phase, the rate of 
baryon number violation is computed to be proportional to exp(— 0(1)4ttv /gT), where v 
is the value of the order parameter, g is the weak coupling constant and T is the tempera- 



ture immediately after the transition [[15| 16| . Thus suppression of the anomalous baryon 



Shaposhnikov has suggested two conceivable ways to enhance the CP violation in the stan- 



dard model at high temperature [45]; the first mechanism, dynamical high temperature sponta- 
neous CP violation, is contradicted by non-perturbative computation [46|, and the second mech- 
anism, reflection of baryon number from expanding bubble walls, according to naive estimates 
cannot provide a large enough asymmetry. 
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number violation after the transition requires a large jump in the Higgs vev v during the 
transition, i.e. 

v(T)/T > 1 @. (1.1) 

Requiring that the baryon asymmetry of the universe be generated during the weak 
transition can give us new information about the CP violating and Higgs sectors of the 
weak interactions, allowing us to rule out some models (such as the minimal standard 
model), and to constrain others. There is still very little known about the origin of CP 
violation or weak symmetry breaking, and any information we can extract from the early 
universe is welcome. 



2. HIGH TEMPERATURE ELECTROWEAK BARYON VIOLATION 

2.1. Baryon Violation and the Anomaly 

Classically conserved global U(l) charges are well known to be potentially violated 
in quantum theories with fermions coupled to gauge fields; this can be seen through the 
one loop computation of the divergence of the corresponding current which leads to the 
anomaly. In the standard model the global baryon and lepton number currents are exactly 
conserved at the classical level. However the anomaly equations give: 

« = « = N f (^-WW - ^XX) (2.1) 

Here Nf is the number of families, is the SU(2) field strength, X^ v the Uy(l) field 
strength, and g and g' are the associated gauge couplings. Note that the difference B — L 
between the two charges is strictly conserved. These equations imply that fluctuations in 
the weak field-strength which have a non-zero dot product of the electric and magnetic 
fields will lead to corresponding fluctuations in baryon and lepton number. The right hand 
side of each of these equations is a total divergence: WW = d^K^ , XX = <9 M /c M where 

K» = f afl (WZ a W$ - \gt ahc W a v WlWl) 

6 (2.2) 

k v = e^X va Xp 

and W, X are the gauge potentials associated with the SU{2) and Uy(1) groups respec- 
tively. Naively this would seem to allow the definition of a new baryon current which 
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is exactly conserved, Jg — Nfg 2 /327r 2 K^ + Nfg' 2 /327t 2 k fl . However is not gauge- 
invariant, and this redefinition is not appropriate |57j . 



If we consider the change in the total baryon number from some initial time zero, 
to some final time r (we will always take space to be a three-sphere, £ 3 ) we can use the 
anomaly to give the change in the baryon number: 

AB = N f [Ncs(r) - N cs (0)} - N f [n cs (r) - n CS (0)} 

Ncs = 32^ / d3xeijk ( W *Wk - l&abcWfWfWS) (2 3) 



/2 p 
9 —r / d 3 X€ ijk XiiX k 



ncS = I a xe ^ij^k 

We have written the result in terms of gauge non-invariant objects Ncs an d ncs, but 
the difference at different times of each of these two objects is gauge invariant. We are 
generically interested in cases where initial and final average values of the gauge field- 
strengths are zero, and we wish to know if AS can be non-zero. In this case ncs is strictly 
zero, since it is proportional to the t/y(l) field strength, which is zero by assumption; 
hence from now on we need only deal with the £77(2) gauge fields and Ncs- To compute 
AS we work in temporal gauge, Wq = 0. Any gauge potential for which the field strength 
vanishes must then be a time independent gauge transformation of W = 0, and so we can 
choose our gauge such that W = Ncs = at t = 0. Then the potential at t = r must be 
a (time-independent) gauge transformation of zero: 

W{t) = hVh' 1 (2.4) 

In terms of h we can write the change in baryon number as 

AB = N f Ncs(T) = N f^2 J d 3 xe ljk Trlhd^hdjh^hdkh- 1 ] (2.5) 

The problem of finding the change in baryon number has been reduced to the classi- 
fication of the possible functions h, which are maps from £ 3 into the gauge group £77(2). 
The parameter Ncs i n this equation is a topological invariant of these maps, known as the 
Chern-Simons number. Thus the baryon number can change if the Chern-Simons number 
can be non-zero. As is well known, the possible values of the Chern-Simons number for 
£77(2) are the integers, and so the baryon number can change by an integral multiple of 
the number of families. 
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2.2. Baryon Violation at Zero Temperature 

Although we have seen that the baryon number can change we have yet to see that 
it does, and we must still compute the rate for baryon violating processes. This requires 
an understanding of the dynamics of the fields as they go from the initial to final field 
configurations. It is this part of the problem which will require non-perturbative methods. 

Non-perturbative effects in gauge theories began to receive attention in the seventies, 
especially in regard to the famous Ua(X) problem, and the problem of baryon violation is 
not too different from the problem of axial charge violation in QCD ||11|| . Since the Uy(l) 
part of the anomaly does not play a significant role, for simplicity we will consider a pure 
SU(2) gauge theory, without the gauged hypercharge. We will also restrict our attention 
to a single family — multiple families will then be easily dealt with. 

If the standard gauge theory of weak interactions does violate baryon number by inte- 
gral units we would expect to find a non-zero amplitude involving three quarks combined 
in a color singlet, plus one electron (or electron neutrino, depending on the electric charge 
of the three quarks). For example we can look at the matrix element (uude); a non-zero 
value for this object would be a sign of baryon violation. In the semi-classical approxi- 
mation this matrix element is given by a path integral dominated by a stationary point, 
the weak instanton. The matrix element above is then evaluated by replacing the fields 
by their values at this stationary point, and multiplying by the determinant given by the 
Gaussian fluctuations about the instanton. So far, everything is identical to the QCD 
case. However, the weak interaction SU(2) symmetry is spontaneously broken, and thus 
the instanton is not an exact stationary point; the Euclidean action coming from the Higgs 
potential can be reduced by shrinking the instanton radius to zero. However, these config- 
urations are approximate stationary points, and should still be the dominant contribution 
to the functional integral. The integral over instanton sizes (including the action from the 
Higgs potential) can be performed to yield a non-zero matrix element which is exponen- 
tially small in the weak fine structure constant a w as is characteristic of any semi-classical 
process: {uude) oc e~ 2ir / aw . The smallness of a w implies that such a matrix element is 
totally unobservable — baryon violation in the form of low energy scattering experiments 
involving small numbers of particles is too small to be observed from this source. 

The generalization to three families is, as promised, trivial. Since each left handed 
doublet has a zero mode in the instanton background, the only non-vanishing baryon 
violating matrix elements must involve at least one fermion from each doublet. Thus the 
change in baryon number during any anomalous event is three times what it would be for 
only one family. 
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2.3. Baryon Violation at non-Zero Temperature 

The semi-classical picture painted in the previous section can be thought of as a 
quantum tunneling: the instanton represents a barrier penetration, and the exponential 
in the coupling is similar to the conventional WKB factor in quantum mechanics. This 
picture can be made precise, by using a canonical Hamiltonian formulation of the gauge 
theory. The barrier penetration factor is then an integral along a path in configuration 
space of the usual WKB functional, which depends on the potential energy of the theory. 

We can get an idea of how this "potential" varies over configuration space by finding 
its minima. We will first consider only the potential energy of the gauge field. One 
minimum is then easy to find: = and no fermions, which we will call Qq. We can 
conventionally take the value of this potential at fio to be zero. Other potential minima 
have gauge fields which are of the form = hd^h -1 . At first sight these are not new 
minima — they are merely gauge transformations of Qq. However if we consider a path 
in configuration space which begins at Qq and ends with = hd^h -1 , we know from 
the anomaly that the final configuration can have non-zero baryon and lepton number, 
both given by the Chern-Simons number of the gauge transformation h. Assuming that h 
does have non-zero Chern-Simons number this configuration is then also a minimum of our 
potential with zero potential energy, but since it has non-zero baryon and lepton number 
is a genuinely different point in configuration space. Consequently the potential has an 
infinite number of minima, which we can label as fi n , where n is the change in Chern- 
Simons number relative to Qq (a gauge invariant quantity); note that this is also just the 
baryon number per family. All that remains is to find the potential at all other points of 
the infinite-dimensional configuration space — an impossible task. However we can plot a 
one-dimensional path in this space which goes through all of these minima, and we choose 
the path for which the height of the barrier between minima is as small as possible. The 
resulting one-dimensional potential is given in fig. [l[ 



The field configuration at the top of the barrier is called the "sphaleron" |13[] , which 
has potential energy: 

E sp = ^B(X/g 2 ) (2.6) 

where B is a constant requiring numerical evaluation; for the standard model with a single 
Higgs doublet this parameter ranges between 1.5 < B < 2.7 for X/g 2 varying between 
and oo (A is the Higgs self coupling). 



Although quantum tunneling through this barrier is irrelevant due to the small cou- 
pling, the finite (although large) barrier height can have important consequences for baryon 
violation. So far we have assumed that the initial and final configurations in our baryon 
violating transitions have no gauge field, but only an assortment of low energy fermions 
with baryon number differing by one; this is the case for the S-matrix element (uude) . But 
we might also consider processes that have more complicated initial and final configura- 
tions, with a total energy greater than the barrier height. In this case quantum tunneling 
is not necessary, as we have enough energy present in the initial configuration to cross the 
barrier. This is the essence of the finite temperature baryon violation argued by Kuzmin 



et al. JL5! — at temperatures near the sphaleron energy initial configurations with energies 
above the barrier will be likely members of the thermal ensemble, this likelihood being 
given by a Boltzmann factor, and baryon number changing transitions can occur without 
barrier penetration. At high temperature baryon number changes by diffusion across the 
top of the barrier, with a rate proportional to exp (—E sp /T). 

We can gain some insight into this problem by imagining a one-dimensional particle 
analogue, a pendulum in the gravitational potential of the earth. This pendulum will be 
the analog of the gauge field (or at least one degree of freedom of the gauge field). For 
the baryons, we imagine that each time the pendulum passes through an angle of tc (the 
point where the pendulum points straight up) traveling in the clockwise direction it moves 
a lever which increases a counter by one; when traveling through n in the counterclockwise 
direction the motion of the lever decreases the counter by one. The configuration space of 
this example consists of the angle of the pendulum (between and 2tc) and the value of the 
counter (B), and the classical potential energy of the pendulum is precisely our fig. []]. At 
zero temperature, barrier penetration will occur, allowing transitions among configurations 
with different values of the counter; if the barrier is high, these processes are suppressed, as 
are the tunneling events in our gauge theory. At finite temperature however, the pendulum 
will experience thermal fluctuations, and at sufficiently high temperature will frequently 
cross 9 = 7r, randomly in both directions. Since each time the pendulum crosses tt the 
counter value changes, we see the value of B fluctuate. If we were to start the system with 
the counter localized at some point in configuration space then B will diffuse away from 
this value like yi, following a typical random walk. 

While this picture does show fluctuations in the counter (and by analogy local fluc- 
tuations in baryon number in the gauge theory case), the average value of the counter 
does not change since fluctuations in each direction are equally likely; this is not yet the 
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physics of baryon violation that we are looking for. As emphasized by Kuzmin et al. |15 
the average value of the baryon number can change only if there is some bias which favors 
fluctuations in one direction over those in the other. The potential of fig. [I] does not show 
any such bias; fortunately this is the result of our leaving out an essential part of the 
physics. We have so far plotted only the potential associated with the gauge field, but 
there is of course an energy associated with the baryons as well. For a sufficiently dilute 
system this is just the energy of a free Fermi gas with fixed baryon and lepton number, 
while at high density it is more complicated. In either case its minimum will be at zero 
baryon and lepton number, and will be symmetric about this value. In our analogy we 
can associate a quadratic potential with the value of the counter, with a minimum at zero. 
Thus there is some curvature associated with the periodic potential, which we show in 
the upper curve of fig. [I]. This additional effect provides the bias we are looking for; if 
we start the system localized near some value of the baryon density, it diffuses outwards 
eventually reaching an equilibrium with zero average baryon number. The bias of the rate 
in the direction of increasing baryon number relative to the decreasing direction is given 



by detailed balance — assuming that this bias is small we may write |2C 



^ = 3(r + -r_) = -3^AF (2.7) 

where AF is the free energy difference between neighboring minima, and T a is the rate per 
unit volume for fluctuations between neighboring minima in the absence of bias (i.e. in the 
absence of fermions). We have assumed that the processes which equilibrate all degrees 
of freedom aside from the baryon number (and the corresponding degree of freedom of 
the gauge field) are rapid compared to the diffusion time, and this is why the free energy 
appears rather than the internal energy. 

2.4. Rates 

The first serious attempt to compute the anomalous rate F a used a technique devel- 
oped by Langer [SS| and Affleck []59[ for evaluating the diffusion rate of a system over a 



barrier at finite temperature. The procedure effectively reduces the theory to one dimen- 
sion (the degree of freedom we plotted above) , computes the flux of the system across the 
barrier in one direction weighted with a Boltzmann factor for each possible initial state, 
and finally reintroduces the degrees of freedom transverse to this one mode in a quadratic 
approximation with a thermal occupation. The idea behind this technique is to treat the 
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most important mode, the one which goes over the barrier most easily, as exactly as possi- 
ble, and then use only the small fluctuations around this mode. The resulting computation 
should be valid when the rate of thermal diffusion across the barrier is small, and the tra- 
jectory over the sphaleron barrier is the most important. The result has all the features 
that we have described so far, especially the characteristic Boltzmann factor associated 
with the sphaleron energy [FTo], [TI]] : 

r a = -f(a w T)- 3 M^e~ Esp/T (broken phase) (2.8) 

Relatively simple arguments have been given for this general form including the exponen- 
tial and the Myy dependence |17|, ^(J . The constant 7 cannot be calculated analytically; it 



depends on the ratio of the Higgs self-coupling to the gauge coupling squared and has been 
numerically evaluated in ref. [IS]. From this form we can see that the Boltzmann factor 



is, naively, a significant suppression up to temperatures of hundreds of GeV. Note how- 
ever that the energy of the sphaleron configuration is proportional to the weak symmetry 
breaking — if the weak interactions were unbroken, we could find a path in configuration 
space which changes the baryon number by one unit, with arbitrarily small potential en- 
ergy at each point. Although the energy in this case is arbitrarily small, the integral of the 
WKB factor along this path is still bounded below by 8ir 2 /g 2 , and consequently quantum 
tunneling is suppressed for any value of the symmetry breaking. On the other hand, the 
thermodynamic baryon violation does not care about the tunneling factor, but only about 
the energy of the configuration, and consequently in the symmetric phase we can change 
the baryon number with neither suppression from quantum tunneling, nor suppression 
from a small Boltzmann factor. This means that baryon violation is only suppressed up 
to temperatures where the weak interactions undergo a transition to the symmetric phase, 
which occurs at a temperature around one hundred GeV. Thus we need the rate of baryon 
violation not where the Boltzmann factor is a large suppression and the calculation of [|19| 
is applicable, but also in the symmetric phase where there is no Boltzmann suppression 
factor at all! 

The calculation of the baryon violation rate in the symmetric phase is considerably 
more difficult than in the broken phase at low temperatures. The problem is that infrared 
divergences begin to invalidate a perturbative expansion when the temperature is higher 
than all mass scales of the theory. These divergences are actually cut off by the genera- 
tion of a "magnetic screening length" £m which gives the scale of spatial correlations in 
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the gauge theory at high temperature, but no reliable perturbative scheme involving this 
screening exists |Kj |6l| . A simple scaling argument allows us to estimate £m- At high 
temperature spatial correlations become increasingly classical, and we need only consider 
the classical dynamics of the gauge field. But in this case the Hamiltonian appears only 
in the combination H/T. Classically the gauge coupling is an irrelevant parameter; by 
rescaling the fields we can always bring the coupling out as a factor l/a w in front of the 
Hamiltonian. Thus the only place the coupling appears is in the combination a w T, and 
therefore on dimensional grounds the screening length is £m ~ {ot w T)~ . Similarly we can 
estimate the rate per unit volume of baryon violating transitions as 

T a = k(ql w T) a (symmetric phase) (2.9) 

where k is a dimensionless constant which requires something beyond a simple scaling 
argument to compute. Ambjorn et al. ||18|| have attempted to evaluate this rate numerically, 
and find a value 0.1 < k < 1.0. These computations are difficult, suffering from both 
infrared and ultraviolet divergences; the results are noisy, and the diffusive behavior that 
one expects is not obvious from the data (for example, growth of the fluctuations in the 
Chern-Simons number as \/t cannot be verified). For these reasons the results must be 
taken as only indicative of the true value. 

Recognizing the difficulties of simulating the standard model in 3 + 1 dimensions, 
Grigoriev, Rubakov & Shaposhnikov [^2| have instead simulated the classical equations of 



motion for a 1 + 1 dimensional analog of the standard model, the Abelian Higgs model. This 
theory has an anomalous fermion number, instantons, and sphalerons similar to the 3+1 
dimensional gauge theory, but is both easier to simulate, and has better ultraviolet behavior 
than its four dimensional analog. Grigoriev et al. used a Metropolis algorithm to generate 
an initial field configuration; since the theory is non-linear and contains an infinite number 
of degrees of freedom, its behavior is hoped to be ergodic, and this initial configuration 
should give "thermal" results, with a temperature associated with the statistical ensemble 
from which it is picked. The classical equations of motion were then used to evolve this 
configuration forward in time. The result of this evolution shows striking confirmation of 
the thermal diffusion across the barrier; as fig. shows, the system spends most of its time 
near one (integer) value of the Chern-Simons number, with only small oscillations around 
this value. However occasionally the system makes a rapid transition to a neighboring 
integer value of Ncsi an d then continues to hover around this new value. The rate of 
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these transitions is in good agreement with the 1 + 1 dimensional analog of (2^8), showing 
the Boltzmann factor associated with diffusion over a barrier. In addition, the rate and 
amplitude of the small fluctuations around one of these integral values is measured to 
be in agreement with analytic estimates. Grigoriev et al. looked in detail at the field 
configurations in the neighborhood of these rapid transitions and were able to verify that 
the field indeed "unwinds" by passing over the barrier between neighboring minima of the 
potential as we have discussed. 

2.5. Baryon Violation at the Phase Transition 

There are several issues which become relevant specifically for baryon violation at the 
electroweak phase transition. If this phase transition is first order and proceeds through 
the nucleation of bubbles of true vacuum with non-zero symmetry breaking amidst a sea 
of false vacuum, we will then need to know the rate of baryon violation not only in the 
symmetric phase (given by ( |2.9| )) or the broken phase when exp(— E sp /T) is small (given 
by ( |2.8| )) but also near the interface of the two phases, the bubble wall. In this case we 
need to know how the baryon violation rate makes a transition between these two forms. 
A naive interpolation between these two formulas as is done in ref. |63|] would indicate that 



the region of the bubble wall where baryon violation is significant would be the region 
where 

M\y (T, </>) < 7a w T (condition for rapid baryon violation) (2-10) 

Here M\y{T, (J)) is the value of the W boson mass which is varying through the wall. The 
rate in this region is still not well known but is likely to be similar to the rate outside the 
bubble. These statements are not definitive since the only known explicit computational 
methods used in ref. |19| are not reliable at this point. 



An entirely different problem has been raised by Turok et al. |34, 65]. As the bubble 



walls created at the weak phase transition expand, they liberate energy in the conversion 
of the false vacuum to the true vacuum. This liberated energy may move the system away 
from thermal equilibrium in the broken phase behind the bubble wall. Then, as is suggested 
by some simulations in 1 + 1 dimensions, none of the formulas above are directly relevant, 
since they assume that most of the degrees of freedom of the system are in equilibrium. 
Turok et al. find that significant Chern-Simons fluctuations persist in this broken phase 
region, contrary to the thermal rate predicted from ref. (|2.8| ). They find that the rate 
for these fluctuations actually increases as the temperature is lowered. However Dine has 
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argued qualitatively that this phenomenon does not occur in the real 3 + 1 dimensional 
world |66[]. At the weak transition the available latent heat is insufficient to create a field 
configuration localized near the top of the sphaleron barrier — the latent heat in a sphaleron 
volume is ~ (n\/g 3 )T, while the corresponding height of the barrier is ~ (47r/g)T, and 
for a transition in which eq. ( |1 . 1|) is satisfied A < g 3 . This analysis disagrees with ref. ||64|| 
where the relevant volume has been taken to be the magnetic screening volume (£m) 3 
rather than the smaller sphaleron volume. This interesting possibility will have to be 
pursued by numerical investigations in 3 + 1 dimensions. 



3. THE WEAK PHASE TRANSITION 

3.1. The Nature of the Transition 

The nature of the weak phase transition in the minimal standard model is fairly well 
understood in two different limits. The transition is known to be strongly first order in the 
Coleman- Weinberg scenario [ |67l - f7l~]| , where the Higgs mass is tuned so that at zero tem- 
perature the curvature of the effective potential for the Higgs vanishes at the origin. By 
"effective potential" , we mean the free energy density for a spatially homogeneous Higgs 
field as computed in perturbation theory by the methods of Coleman & Weinberg |57, 
72fl . For values of the Higgs field between two local minima this is not the "true" effec- 



tive potential, which includes spatially inhomogeneous field configurations and is always 
convex []73|, [74| . The Coleman Weinberg effective potential in principle is relevant for com- 
puting properties of the phase transition such as the existence of metastable phases and 



bubble nucleation rates [58, 72]. In the minimal standard model the Coleman Weinberg 
scenario requires a Higgs mass of 10 GeV, which is ruled out; however in models with 



several scalars this scenario is still viable [ 75] , [f6| . As the mass of the Higgs increases, the 



strength of the transition decreases. Perturbative calculations ]45], |68|, [71], |7q - |88[1 indicate 
a first order transition persists as the Higgs mass increases. However in the minimal stan- 
dard model the lower bound on the Higgs mass is currently 60 GeV, comparable to the 
W mass ||89j] . For Higgs masses greater than the gauge boson masses finite temperature 
perturbation theory is unreliable at the critical temperature El], |83|-|88|, |90r|93j] . In the limit 



where the Higgs mass is much heavier than the W and Z masses, it is a reasonable approx- 
imation to neglect the gauge and Yukawa interactions when studying the transition. Then 
the minimal standard model is approximately the linear sigma model for 0(4) symmetry 
breaking. The finite temperature phase transition for the 0(4) model has been studied 
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numerically using lattice field theory [71, 9|, 95] and by imposing an infrared cutoff which 
is removed using the renormalization group |)6], and is clearly second order. Thus for a 
heavy Higgs the transition is at most very weakly first order, as is also verified by lattice 
calculations including the gauge couplings |97|]. Computations with N flavors of scalars in 
the large N limit confirm the result that the transition is second order or weakly first order 

1. It has also been 



for a heavy Higgs, and first order for a sufficiently light Higgs |77], |98|, 
argued using the renormalization group (see ref. ||100|| ) that the transition is first order for 



small Higgs mass and second order for large Higgs mass ||101 



The weak phase transition is under active investigation since demonstrating the 
feasability of electroweak baryogenesis and doing a precise computation of the resulting 
baryon abundance requires an understanding of the details of this transition — the order, 
the critical temperature, the bubble nucleation rate, the velocity and shape of the expand- 
ing walls, and the Higgs vev after the transition. Most of these details depend on the model 
of the weak symmetry breaking sector, about which we currently have very little experi- 
mental information. Requiring a first order transition with a sufficiently large Higgs vev 
after the transition places significant constraints on this sector. For instance, in the mini- 
mal standard model the condition ( |1 . 1|) for avoiding washout of the baryon number after 



the transition places an upper bound on the Higgs mass [55]. The numerical value of this 
bound is still uncertain, but perturbative calculations place it in the experimentally ruled 
out 30-45 GeV range [[45], |8(], |4|, |k| |102| , |103|| . In models with additional scalars there is 



a weaker upper bound on the mass of the lightest scalar [ff6| , [81] , |104H110|| . For instance 
Anderson & Hall showed that simply adding a gauge singlet scalar s with a s 2 \H\ 2 
coupling substantially weakens the Higgs mass bound. Similarly, Myint ||106|| showed that 
in the minimal supersymmetric standard model, the coupling of the top squarks to the 
Higgs doublet weakens the Higgs mass bound slightly. 

The main theoretical tools for studying the transition are finite temperature perturba- 
tion theory and lattice field theory. The lattice studies of the finite temperature transition 
for most of the interesting cases, such as for the minimal standard model including gauge 
boson effects with a relatively light Higgs, are still in a preliminary stage [71, [J7J, |1 1 1| . 
112(1 , and do not contradict improved one loop perturbation theory []77|, |75|, |8"5| - |8"S| , |113| - 



115[| . Computing the effective potential for the order parameter </> gives the order of the 



transition — if at the critical temperature there are two degenerate minima the transition 
is first order. In fig. |3| we plot the effective potential at various temperatures for a typical 
example of a first order phase transition. 
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The perturbative one loop result for the effective potential V is ]67j, |68|, 

1 



V{<f>,T) =Ki(0) + 
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& 2 d£; 



(2tt 
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exp 



(3.1) 



where Mi and are the mass and spin respectively of the z'th particle, the sum on % includes 
all degrees of freedom including spin, V c \ is the classical potential and the polynomial 
in (j) is determined by the choice of renormalization condition. Physically this result is 
interpreted as the ground state energy plus the free energy of a gas of noninteracting 
particles at finite temperature. We therefore do not expect eq. ( |3.1| ) to be accurate at 
arbitarily high temperatures, where the density of thermal particles is too high to neglect 
interactions. Examination of higher order corrections for a theory with gauge interactions 
and scalar self coupling A shows that the naive finite temperature perturbation theory at 
high temperature is an expansion in powers of XT 2 /M 2 ((p) and gfT 2 /M 2 ((p), where gi is 
the coupling of the i'th gauge particle with mass Mj to </>. Clearly this expansion is not 
valid at high temperatures, or for any values of <p for which there are massless interacting 
particles. However this one loop result can be improved at high temperatures, by summing 
up the dominant high temperature corrections to all orders. The leading corrections, which 
grow as T 2 , come from subgraphs with superficial degree of divergence 2 — i.e. corrections 
to the boson propagators [[77], [7JJ . How best to improve perturbation theory is currently a 



subject of active investigation |83| - |86| , |113| - |119|| . Here we will summarize the basic idea of 
modifying the boson propagators to sum a class of diagrams which give the leading high 



temperature corrections. For T much larger than Mj the integral in eq. (|3.1|) becomes 
approximately 

boson M 2 (0)T 2 



for a boson and 



fermion 
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for a fermion, and gives an effective temperature dependent contribution to the scalar mass 
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Including this mass in the scalar propagator is equivalent at one loop to the summation 
of all the leading high temperature contributions of the "daisy" graphs such as those in 
fig. |] [[77, 75] . I n perturbation theory there is also a temperature dependent contribution 
to the effective mass of the longitudinal gauge bosons — the usual Debye mass of order 
gT. The transverse components of the gauge bosons are believed to get a non perturbative 
contribution to their masses of order a w T, but remain massless to all orders in perturbation 
theory. A systematic way of including the leading 0(T 2 ) high temperature contribution 
to the 4> propagator []87], |115|1 is to add and subtract from the Lagrangian a term 

C - C + y AM| eff - ^AM| eff , (3.5) 

where the first term is treated as part of the unperturbed Lagrangian and included in the 
scalar propagator while the second term is treated in perturbation theory as an interaction 
which cancels the leading high temperature pieces of thermal subloops. A similar but more 
complicated method can be used to sum the leading high temperature contribution to the 
gauge boson propagators (see ref. |87[]). Improved perturbation theory then becomes an 
expansion in powers of AT '/ 'M </,((/)) ff and g 2 T/Mi(<fi) eS , up to logarithms. Because the 
effective thermal masses of the particles grow linearly with T (for large T) the improved 
expansion does not necessarily break down at high temperature. It always breaks down 
near a second order phase transition however, since M^(0) eff goes to zero as the critical 
temperature is approached, and is also unreliable in the vicinity of a very weakly first 
order transition. For a sufficiently strongly first order transition, it is possible for the 
expansion parameters to remain small in both minima of the effective potential at the 
critical temperature. 

For the minimal standard model with a Higgs mass in the (30-80) GeV range, a decent 
approximation to the improved one loop effective potential is to compute eq. ( |3.1[ ) using 
the temperature dependent masses for the bosons and a high temperature approximation 
for the integrals, the result is [B4f : 



M A^ 4 _ (2„4 + m|)r + (2m^ + raj + 2ra^ - Wf + covst , (3.6) 
4 67ti>q 8vq 

where 



_2v 2 [ml - 3(2m^ + m % - Amt)/(^ 2 vD] 
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Note the presence of the \(f>\ 3 term, which means that at the critical temperature there will 
be two coexisting minima, hence a first order phase transition |6?|] . For the range of Higgs 
masses considered, this critical temperature T c will be just above T . Is this calculation 
reliable? At the critical temperature, the effective potential has the form 



[2m 3 w + m%)T t 1 2 



3ttv$\ Tc 



(3.f 



Note that the condition (|1.1| ) requires 



A Tc < 9 3 ■ (3.9) 
The effective scalar mass in each of the two minima is 



, , 1 (2m 3 w + m%)T c 
M ^ ~ V 2A7 3^1 • ^ 

In either minimum, additional loops involving the scalar self coupling will be suppressed 
by 



\ T T C V 2A ^ 37r ^ m H 



(3.11) 



M^ eff (2m 3 w + m%) (2m 3 w + m%) ' 

The gauge bosons in the broken phase and the longitudinal gauge bosons in the symmetric 
phase are heavier than the scalars and will give rise to corrections proportional to pow- 
ers of g. However the transverse gauge boson masses in the symmetric phase vanish in 
perturbation theory; perturbative corrections involving these states become unsuppressed 
for < gT and are infrared divergent, making a perturbative calculation of the effective 
potential at the origin impossible |5D[ [JT|. If we assume that these modes receive a 
nonperturbative mass of order a w T (related to the magnetic screening length), the uncer- 
tainty in the value of the effective potential at the origin is 0(g 6 T 4 ) due to nonperturbative 
effects []60| , |6T| . Nevertheless it is still possible to establish the order of the transition in 
perturbation theory. If At c is small (light Higgs mass) when compared with g 2 then ex- 
amination of eq. ( |3.8| ) shows that the barrier between the two minima occurs at relatively 
large (p and the height of the barrier is larger than the g G T 4 uncertainty in the value of the 
effective potential in the symmetric phase. (The calculation of the effective potential also 
breaks down near the inflection points, but a lower bound on the barrier height can be in- 
ferred from the value of the potential in the region between the minimum and the inflection 



18 



point where the perturbation expansion parameters are small.) Thus while the nonpertur- 
bative effects will induce some uncertainty into the value of the critical temperature, they 
will not affect the result that the transition is first order provided that 

gT c « 1 f 3 . z) , (3.12) 

which implies that m# < raw- In the minimal standard model with a Higgs mass in the 
experimentally allowed range eq. ( |3.12| ) is at best marginally satisfied, and so perturbative 
calculations should not be terribly accurate. More sophisticated calculations validate this 
conclusion — for instance Arnold & Espinosa |87| and Bagnasco & Dine |88|] find fairly 



large (40%) two loop corrections to the improved one loop result for the value of the Higgs 
vev at the critical temperature in the minimal standard model with a 60 GeV Higgs. 
Clearly better methods are needed for quantitative results, such as nonperturbative lattice 
calculations. 

3.2. Dynamics of the Transition in the Early Universe 

The conventional picture of how a first order cosmological transition proceeds is that 
as the universe expands it cools in the symmetric phase until it reaches the critical temper- 
ature T c at which time the broken phase becomes equally favored. At T ~ T c bubbles of 
the broken phase nucleate, but surface tension effects cause these to immediately shrink, 
so that the universe remains in an approximately homogeneous state. For T somewhat 
below T c , bubbles of broken phase nucleate which are sufficiently large so that the volume 
pressure from the lower free energy inside the bubble can overcome the surface tension. 
These critical bubbles grow until the universe is completely converted to the broken phase. 
The time scale for the formation of a critical bubble is proportional to exp(i ? c /T), where 
F c is its free energy. At T = T c , F c = 00 and so the universe always supercools, resulting 
in a departure from thermal equilibrium ||58|| . A detailed study of the bubble growth in 



the minimal standard model is given in ref. ||120 



This picture was recently challenged by several authors who claimed that thermal 
fluctuations in the order parameter could maintain thermal equilibrium even during a 
first order transition [|9~0"l - |92"ll , which would make baryogenesis impossible. For instance 



subcritical bubbles of the broken phase could appear above the critical temperature, and 
during the transition the universe could be an emulsion of both phases with the fraction of 
broken phase gradually increasing. Such large thermal fluctations would also destroy the 
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validity of perturbation theory, since perturbative calculations assume that the deviations 
from a homogeneous field configuration are small ]84], |93|, 121|| . Fortunately, for those 



models in which eq. (|1.1|) is satisfied such thermal fluctuations are suppressed, perturbative 
expansion parameters are less than 1 , and the conventional picture is probably correct [|S4| , 

n nn HI- 

Another conventional assumption is that the bubbles are nearly spherical and that, 
once they have grown to macrosopic size, it is a good approximation to take the walls to 
be planar. It has been suggested [|123|1 that smooth walls are unstable towards developing 
wrinkles, which would greatly complicate bubble evolution, but a recent computation | ]124| 
indicates that for the weak transition smooth walls are stable and that surface tension will 
keep the expanding bubbles spherical until they collide. 



3. 3. Propagation and Shape of the Bubbles 

In the conventional picture of the transition local departure from thermal equilib- 
rium, and baryogenesis, primarily occurs in the vicinity of the expanding bubble walls. 
A quantitative calculation of the baryon number of the universe as described in the next 
section requires the temperature at which bubbles nucleate, the wall profile, and the wall 
velocity. If the effective action is known, computing the nucleation temperature and the 



wall profile is a straightforward exercise [|12| , |5q , p9|, |125| . The nucleation temperature 
is the temperature where the rate per unit volume of critical bubble nucleation equals 
the age of the universe multiplied by the volume inside the horizon. The critical bubble 
nucleation rate depends on F c , which can be computed by finding a spherically symmetric 
field configuration <p(r) which is a stationary point of the effective action and has </> = 
for large r, and </> near the minimum of the effective potential for small r. The free energy 
for such a configuration is approximately (neglecting corrections to the derivative terms) 



S[<t>(r) 



r dr 



L'\2 



(3.13) 



and a stationary point can be obtained by finding a nontrivial solution of 



r 2 dr dcj) 



(3.14) 



satisfying the boundary conditions 



cj)(oo) = o, cj)(oy = o . 



(3.15) 
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This must be done numerically, although an approximate analytic formula is given in 



ref. [S4|. In the minimal standard model with a Higgs mass in the 30-80 GeV range the 



result is that the nucleation temperature is nearly the critical temperature [|63| , |ST| , |8~2 



103|| . The wall profile for the boundary between the two coexisting phases at the critical 
temperature is the solution to the differential equation ( p. 14] ) at T = T c in the large r limit, 
subject to the boundary conditions 4>(oc) = and <f>(— oo) = v, where v is the vev in the 
broken phase. This solution should also give a good estimate for the profile of the wall of 
an expanding bubble. (The exact wall shape of an expanding bubble in a thermal medium 
is still unknown.) If the effective potential at the critical temperature is (A/4)(0 — v) 2 <p 2 , 
the analytic solution for a wall centered at r = is 



v 
2 



1 + tanh ( — r v 



(3.16) 



The width of the wall 5 W = (y/S/X )/ v is quite model dependent, which is unfortunate 
since, as we show in §4, this is an important parameter for calculating the baryon num- 
ber produced. In the minimal standard model, if one takes the Higgs mass to have the 
experimentally ruled out value 35 GeV so that v/T fa 1, the wall widthi is 5 W ~ 24/T. In 
general 5 W oc vf yh, where h is the barrier height between the two minima of the effective 
potential. In a model with a strongly first order transition, such as a version of the singlet 



majoron model with the potential tuned to have a nearly flat direction |76], or a super- 
symmetric model with a gauge singlet superfield [p. 101 , the wall can be thin, 8 W = 0(1/T). 
These calculations require knowing the effective potential in the region between the two 
minima, and perturbation theory breaks down at the inflection points. In between the in- 
flection points the perturbative potential has an imaginary part, arising from the classical 
instability of the homogeneous field configuration which we have assumed dominates the 



perturbative expansion |7^|. However the lifetime of such a homogeneous configuration is, 



in most cases, long enough that it is meaningful to compute its free energy, which is ap- 
proximately given by the real part of the effective potential |7j| . Any smooth continuation 



of the effective potential from the regions where perturbation theory is sensible will not 
deviate too far from the perturbative value, and for a sufficiently strongly first order tran- 
sition in a weakly coupled theory the perturbative calculations should give a reasonable 
estimate for nucleation rates and the wall width <$„,. 



3 This estimate differs from ref. |84| because of a factor of two difference in the definition of 
the wall width. 
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The calculation of the wall velocity v w is more complicated, since it involves nonequi- 
librium interactions of the thermal plasma with the wall. The higher pressure inside the 
bubble will tend to accelerate the wall, and the scattering of thermal particles from the 
wall will dissipate energy and slow the wall down [125 ] . A plausible assumption is that the 
wall will reach a terminal velocity where the different forces balance. This assumption was 



challenged by Turok [|126|, who pointed out that if the temperature is spatially constant 



and the particle distributions are always the local thermal distributions for a given value 
of <fi, then the force on the wall is independent of velocity and continuously accelerates the 
wall. This is of course true, since the wall can only slow down by dissipating energy, which 
will cause a local departure from equilibrium. Thus if the wall is extremely thick when 
compared with particle mean free paths, so that the particle distributions approximately 
equilibrate with the changing value of (p, v w will be ultrarelativistic. For a wall which is 
not infinitely thick, particle interactions with the moving wall will affect their distributions 
in a velocity dependent way. For instance weak gauge bosons and top quarks, which are 
massless in the symmetric phase and have large masses in the interior of the bubble, will 
tend to reflect from the wall. The reflected particles absorb some of the wall's momentum 
and the distribution of reflected particles is nonthermal. Several estimates of v w , which 
include the effects of particle interactions with the wall and rethermalization in various ap- 
proximations, indicate that v w ~ 0.1c for a wall of width 5 W ~ 1/T and is at most mildly 
relativistic for 5 W ~ (10 — 30)/T [j63] , |84] , |127|| . These calculations assume the minimal 
standard model phase transition parameters with only 5 W allowed to vary; in general v w 
is model dependent. 



4. CALCULATING THE BARYON ASYMMETRY 

The master equation for weak scale baryogenesis is eq. fl2.7|) , which may be rewritten 



as 

dt T dB T K ' J 

where we have replaced AF = 3dF/dB, the derivative computed at fixed B — L, remem- 



bering that AB = 3 in an anomalous event. The rate T a is given in eqs. (|2.8|) , (|2.9|) . We 
have given above the reasons to suspect that two of the three necessary ingredients for 
baryogenesis may exist at the weak phase transition: baryon violation and departure from 
thermal equilibrium. There remain the requirements of C and CP violation, as well as 
some concrete mechanism that explains how a nonzero \i B can arise dynamically during 
the weak phase transition. 
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4-1. CP violation 

We wish to explain the observed baryon to entropy ratio n B /s ~ (0.6 — 1) x 10 -10 
0]. From the expression (|2.9|) for the anomalous rate T a we expect a factor of ~ 10~ 6 
in our final answer for n B /s, and since s = (2n 2 g*/4:5)T 3 with g* ~ 10 2 at the weak 
scale, one can only tolerate an additional suppression of order 10 -2 from the combined 
effects of CP violating angles Sep, and various dynamical effects. This guarantees that 
CP violating effects must be large at weak scale energies, further constraining possible 
models and giving reason to hope that new CP violating effects may be experimentally 
accessible beyond those seen in the kaon system. 

It is natural to wonder whether the requisite CP violation at the weak scale could be 
of the Kobayashi-Maskawa (KM) form. At high temperature one may expect that that 
a dimensionless measure of CP violation 5c p in the standard model could only depend 
on a reparametrization invariant combination of the dimensionless Yukawa couplings and 
mixing angles. In §1 we argued that Sep — 10 -20 , a suppression factor apparently ruling 
out KM CP violation as playing any role in weak scale baryogenesis. In the kaon system, 
the CP violating parameter e is seventeen orders of magnitude larger than this estimate 
because the CP conserving contribution to K° — K° mixing is so small. It is conceivable 
that dimensionless numbers such as m p i/T or S W T could play a role in enhancing Sep the 
necessary eighteen orders of magnitude during baryogenesis. Attempts in this direction 
have been made in ref. [^3J , but the existence of these enhancements remains unconvincing. 
It is most likely that any model of baryogenesis requires new sources of CP violation beyond 
the standard model, and that electroweak baryogenesis requires the new CP violation to 
be at the weak scale. 

Incorporating additional CP violation into the standard model requires additional 
matter fields, and there are a number of relatively simple extensions one can consider. 
Three examples that have been analyzed in the literature and are discussed below are 
the singlet Majoron model with complex phases in the neutrino mass matrix; the non- 
supersymmetric two Higgs doublet model with a CP violating relative phase between the 
doublets; and the minimal supersymmetric standard model (MSSM) with a CP violating 
phase in the gaugino masses. In each of these models 5c p can be 10 -2 or larger, not being 
suppressed by small couplings as in the KM model, and without being in conflict with 
experimental limits on the electric dipole moments (EDM) of the electron or neutron. 
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4-2. Timescales 

Having assembled the ingredients is not enough to make a universe: we still require a 
recipe. The recipe involves the three distinct time scales governing the relevant reactions. 
These are the thermalization time r T which characterizes how fast particles in the cosmic 
plasma equilibrate; the Higgs time scale r H which governs the departure from equilibrium, 
given by r H ~ (f>/<p as measured by a comoving observer while the expanding bubble wall 
passes through her; and the sphaleron time r sp which governs the rate of baryon violation 
in the symmetric phaseH 

Since baryon violation in the broken phase within the bubbles must be very slow to 
prevent re-equilibration of the baryon number, we will assume here that baryon number is 
exactly conserved inside the bubbles. Outside the bubbles in the symmetric phase, r sp at 
the transition temperature T is estimated as 

r" 1 ~ a* T = 1.3 x 10" 6 T . (4.2) 

In contrast, the thermalization rate due to weak and strong interactions, defined for the 
different particles as the inverse mean free path £ T , is much faster, not having the 
suppression. It is estimated from strong or weak Coulomb scattering cross sections: 

r" 1 = r 1 ~ |°- 25T q (4 3) 

These rates correspond to "fast" reactions; some standard processes are much slower, of 
course, such as chemical equilibration between the first and third families, or between the 
two chiralities of a light fermion. 

The Higgs time scale is less easily determined, and is quite model dependent, roughly 
given by 

r" 1 = ~ ^ ~ (0.01 - 1.0)T . (4.4) 

o w 

We see that baryon violation is always out of equilibrium near the wall, since r sp ^> r H . 
However, other particle interactions may or may not be able to equilibrate near the bubble 
wall, depending on the relative size of r H and r T . This gives rise to two distinct regimes: 

1. The adiabatic regime: r T Ct h . "Fast" interactions maintain thermal equilibrium as 
the bubble wall passes by and the value of the Higgs field changes. This allows one to 

In principle there is a fourth time scale, the age of the universe T/f = H~ 1 , but at the weak 
epoch it is many orders of magnitude larger than the others and not relevant for microphysics. 
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to describe the plasma within the bubble wall in terms of equilibrium thermodynamics 
with quasistatic chemical potentials for quantities that equilibrate slowly compared 
to t h . In this regime baryogenesis occurs within the bubble wall itself. 

2. The nonadiabatic regime: r T ^> r H . The wall is thin compared to the mean free path 
and particles reflect off the oncoming wall with calculable CP violating reflection 
coefficients. Baryogenesis occurs in an extended region preceding the phase boundary 
and is enhanced relative to the adiabatic scenario. 

Which regime one is in depends on the bubble wall thickness and velocity: for thick or slow 
walls it is the adiabatic regime that is relevant; for thin or fast walls it is the nonadiabatic. 

4-3. The adiabatic "thick wall" regime: spontaneous baryogenesis 

In the adiabatic limit it should be valid to treat the plasma in the bubble wall as being 
in quasi-static thermal equilibrium with a classical, time dependent field. The plasma will 
not be in chemical equilibrium, however, as some reactions, such as baryon violation, are 
slower than r~ . This deviation from chemical equilibrium may be treated by introducing 
chemical potentials for the slowly varying quantities. 

Baryogenesis arising from slowly varying classical Higgs fields during a phase transition 
is called "spontaneous baryogenesis" ; it was first discussed in refs. ||128| , |1J in the context 
of second order phase transitions, for which the assumption of adiabaticity is always valid 
(except in the extremely early universe). To see how time dependent Higgs fields can 
drive baryon production consider a toy model with a conserved £^(1)^ symmetry carried 
by fermions if) {B = 1), x (B = —1), a scalar <fi (B = 2), and the Yukawa interaction ifxpx- 
During a phase transition in which <p takes on the classical value 

4>(t) = p(t)e"« , (4.5) 

this leads to baryon production, even for (p small compared to the fermion masses. To see 
this, redefine the fermion fields by the phase e z9 ^ B / 2 - : this removes the phase from the 
Yukawa interaction, but leads to a new interaction from the fermion kinetic term: 

Ck.e. -> Ck.e. + ^(^7> - X7°X) = Ck.e. + \6n B ■ (4.6) 

Thus 9 acts like a chemical potential for baryon number. (We use the word "charge poten- 
tial" for 9 instead of chemical potential, since the effect is dynamical and not arising from a 
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constraint.) The interaction (|4.6| ) splits the energy levels of baryons versus antibaryons so 
that the free energy is minimized at nonzero baryon number. The interaction (|4.6|) violates 
C and CP and will spatially average to zero except in theories with explicit violation of 
these symmetries. 

If B symmetry is violated in some other sector of the theory {e.g., through the elec- 
troweak anomaly), then the interaction ( f4.6j ) will cause B to try to equilibrate to the value 
n B = 0{9T 2 ). If the B violating interactions are rapid compared to 9/9 then this equi- 
librium value will be attained. In the real case we wish to study, however, the rate per 
unit volume of anomalous fluctuations T a { see eq. ( |2.9|) ) is too slow to equilibrate due to 
the factor of a%. In this case one simply integrates h B in eq. (|4.1| ) over the appropriate 
timescale. It should be emphasized that it is not necessary for the charge potential to 
couple to the baryon current as in the toy example Q4.6|) in order for there to be a nonzero 
value for h B . As long as the charge X is not orthogonal to B, an interaction coupling of 
the form 9j will give rise to a nonzero value for \x B of the form 

fi B =M9, (4.7) 

where M is a calculable, model dependent constant. 

After the bubble wall has passed by a given point in space, the charge potential 9 
returns to zero. If baryon violation is still rapid compared to the cooling rate of the uni- 
verse, the baryon number produced will re-equilibrate to B = 0. However, the anomalous 
fluctuation rate T a goes rapidly to zero as the Higgs vev turns on, so that the baryon 
number drops out of chemical equilibrium and remains to the present epoch. By means 
of the inequality ( |2.1U| ) we can make a crude estimate of this "cut-off" value of the Higgs 
vev: 

(0} co ^ ^ = 0.7T . (4.8) 



To find the total baryon number produced during the phase transition, eqs. (|4.1|) , (|4.7f) 
may be integrated across the bubble wall: 

Q/V f tco . QA/T 

U ' dt9T a {(j){t)) ~ -— , (4.9) 



where we have estimated the integral by treating T a as a step function equal to its sym- 
metric phase value for {(p) < {4>) co and vanishing further inside the bubble. The quantity 
A9 is then the change in 9 from the symmetric phase to the point in the bubble wall where 
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baryon violation effectively shuts off; it is a homogeneous function of the CP violating 
parameter 5c p of the theory in question. Furthermore, AO is homogeneous in (0} co and 
could be very small if (4>) co proves to be as small as suggested in refs. |5T| , |103| . 

Given that the entropy density is s = (27r 2 g*/45)T 3 and T a = /t(a ffi T) 4 , one finds for 
the final baryon to entropy ratio 

?f = -k (J^J (j?) x (3 x lO" 7 ) (adiabatic) , (4.10) 

to be compared to the observed value n B j ' s ~ 10 -10 . In the above expression, the quantity 
M defined in eq. (|4.7|) is typically 0.1 — 1.0; k is the prefactor appearing in equation ( |2~9| ) 
and is thought to be 0(1); (100/g*) is also close to unity at the weak phase transition. 
The reason why this result is a formula for the baryon density of the entire universe, even 
though baryon production is basically a surface effect only occuring within the bubble 
wall, is that the entire universe is swept out by the bubbles during the course of the phase 



transition 48 



We now list the steps that must be followed when analyzing weak scale baryogenesis 
models in the adiabatic limit, and then proceed to give several examples. 

(i) Identify a model that has a first order phase transition with large timescale r H and 
classical Higgs field evolution during the phase transition such that there is a CP 
violating, time dependent phase 9{t) in the Yukawa interactions. 

(ii) Determine A9 from the Higgs equation of motion in order to use eq. ( |4.10| ) to compute 
the final baryon number. Using the generous estimate ( |4.8j ) for (0) co , A9 will typically 
be 0(1) x 5c p- Additional suppression is possible if baryon violation shuts off closer 
to the outer surface of the bubble and (4>) co proves to be much smaller than assumed. 

(iii) Rotate the fermions in the theory to remove the time dependent phase from the 
Yukawa interactions into a derivative coupling of the form 9j . The current X is 
ambiguous, since the rotation is not unique — one can always rotate so that X — ► 
(X + Q) , where Q is a classically conserved charge. It is simplest to choose X so that 
the rotation is anomaly-free, so one need not treat 6{t)WW couplings. 

(iv) Compute the dynamically generated \i B . To do so one must introduce chemical po- 
tentials for each conserved and appoximately conserved quantum number. All the 

27 



particle densities may be determined in terms of these chemical potentials and the 
charge potential 9; for particle species i the density rii is given by 



ki 



Xi0 + X>a<7? 



rp2 

~6~ 



(4.11) 



where ki is a statistical factor accounting for effective degrees of freedom and differs 
for fermions and bosons, Xi is the X charge of species i, and the sum is over all charges 
out of chemical equilibrium (including B). The values of the chemical potentials (j, a 
are fixed by requiring that the primordial plasma not carry any quantum numbers. 
One can then determine \i B and the factor Af in eq. Ql7| ). 

(v) Use the expression ( f4.10|) to give the final baryon to entropy ratio. This formula is 
only expected to be valid deep in the adiabatic thick-wall regime, however. If the 
equilibration time for particle distributions is long compared with r H there may be 
additional supression factors. 

We now present two examples of standard model extensions and show how spontaneous 
baryogenesis proceeds in them. The first is the two Higgs doublet model, and the second 
is the minimal supersymmetric standard model (MSSM). 

The two Higgs model. This model was first suggested by McLerran as a model for weak 
scale baryogenesis in ref. [ 129|| and analyzed by Turok and Zadrozny in ref. [49]; the 
treatment here will follow our subsequent work f5^] . The model is characterized by the 
Lagrangian 



fermions 



/If! | 



V((f>i, 02 ) + £yuk + £gauge 



i=l 



where the scalar potential is given by [130 



V{fa, 0a) = Xi((p\(pi - v\f + \ 2 {4<h ~ 4) 2 + A 3 [(0l0i - ^i 2 ) + (402 - vl)} 2 
+ A 4 [(0i0i)(0|0 2 ) - (0 t 1 2 )(0j0 1 )] + A 5 [Re(0l0 2 ) - v lV2 cos^] 2 
+ A 6 [lm(0j0 2 ) - wi^ 2 sin£] 2 . 



(4.12) 



(4.13) 



The Yukawa interactions in £ yu k couple the 0i field to up-type quarks; down-type quarks 
and charged leptons may be coupled exclusively to either 0j or 2 — either choice suppresses 
flavor changing neutral currents and is protected by a discrete symmetry that is softly 
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broken. For convenience, we will assume that all fermion masses arise due to interactions 
with (pi ; the other option works equally well, since only the top quark Yukawa interaction 
proves to be relevant. In either case we can adjust the phases of the Higgs doublets so that, 
in the fermion mass eigenstate basis, the top quark Yukawa interaction y t is real and <\>\ 
develops a real expectation value at zero temperature. In this basis, aside from the usual 
KM phase, there is an additional CP violating phase £ in the scalar potential, so that the 
zero temperature vev of the fields have a relative phase e 1 ^ between them. During the weak 
phase transition the neutral components of the two doublets will acquire time dependent 
classical values with 4>\{t) = p(t)e l9 ^ (in unitary gauge where the linear combination of 
the phases of (jy{ 2 eaten by the Z has been rotated away) . The change in 9 in going from 
the beginning of the phase transition to T = is ||52|| 

A6 = o(t- arctan ^ tan£ V (4.14) 

A9 will be smaller in going from the symmetric side of the domain wall to the point where 
the scalar vev equals (4>) co and baryon violation turns off, since A9 scales like ((</>) co /v) 2 , 
arising from the change in relative importance between the <j) 2 and cf> terms in the potential 

(EH)- 

Next we remove 9 from the Yukawa interactions by performing a time dependent 
rotation of the fermion fields in the Lagrangian. Since we do not wish to complicate our 
task and induce any coupling of 9 to gauge fields, we eliminate 9 by means of an anomaly 
free fermion rotation; thus we rotate the fermions by an amount proportional to their 
hypercharge. We find that 9 now has derivative couplings to twice the fermionic part of 
the hypercharge current (the "X" current in this model) 

£ = %e ~ V(<l>i, + VtQLPitR + VblLPi^R + ■■■ + h.c. 

_ _ _ _ (4.15) 

+ 2<9 M [\q L ^q L + \U R -fU R - \D R -fD R - \l Ll n L - Er^Er] , 

where there is an implicit sum over the three families in the current coupling to 8^9, and p\ 
is the radial component of (j>\. The quantity 9 which is nonzero during the phase transition 
is then a charge potential which splits fermion-antifermion energy levels proportional to 
their hypercharge inside the domain walls. 

We now must calculate the thermodynamic force \x B = MO responsible for the pro- 
duction of baryon number. To this end we introduce chemical potentials for each quantum 
number that equilibrates slowly compared to r H . The relevant ones for this model are 
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B — L, Q, and B% — \{B\ + B2) where Bi is the baryon number of the i th family. (There 
are other conserved or slowly varying charges, such as B1—B2, but their chemical potentials 
turn out to be irrelevant.) Prom eq. ( [4.1 1| ) we have determined 



Af 



6 + 71 

3 \ 25 + An 



(4.16) 



where n is the number of physical charged scalars which are light compared to the tem- 
perature T; for n = 0, Af = 0.3. Comparing eq. ( |110|) with G4.14j ) and (|4.16|) , we find 
that a value of £ ~ O(10~ 2 — 10 -3 ) for the CP violating parameter in the theory yields 
an acceptable value for the final baryon to entropy ratio produced in the phase transition. 

Our treatment is justified as long as fermionic hypercharge has time to adjust itself 
through the top quark mass or Yukawa interaction with thermal Higgs particles. This 
adjustment can occur via scattering (e.g., t + g — > t + </>i), scalar bremstrahlung (e.g., 
t+q — > t+q+(pi) or coherent top mass oscillations within the domain wall. Comparing with 



photon bremstrahlung ||131|| we arrive at a crude estimate giving a fermionic hypercharge 
equilibration rate of Ty 1 — t" 1 x (X t /n) 2 ~ 0.03T, where r" 1 is the strong interaction 
thermalization rate (fOf). Thus the above analysis is expected to be valid for r" 1 at the 
low end of the range ( fOj) . For less adiabatic phase transitions there may be an additional 
suppression of our result proportional to (\ t /7r) 2 due to the inability of the plasma to reach 
complete chemical equilibrium in the time available. 

A different treatment of the two Higgs model is found in refs. f£3], [5t|, ||64|| . They 
find a contribution to the effective action in an expansion of rrit/T of the form 



AS 



<(4 ( ^ 

2 \TTVi 



ts ld x 



(Vrfl^Vofa + h.c. 



ijk rpa 



(4.17) 



which is argued to bias anomalous baryon violating transitions when the Higgs field <pi 
takes on the classical value ( [4.5| ) ; the computation is shown in detail in ref. |54j , where an 
estimate of n B /s ~ 10 -10 x (CP violating angles) is given. This computation is subleading 
in the adiabatic limit, as it is proportional to (At/7r) 2 . Furthermore the effect of this 
operator is more strongly affected by the value of (<fi) co where baryon violation cuts off, as 
its effects scale like (<ft) co [p4" l, while the factor of A# in formula (|4.10| ) scales like (4>) co - 
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Minimal Supersymmetry. One of the most attractive extensions of the standard model 
is the minimal supersymmetric standard model (MSSM) ||13Q|] . It possesses two Higgs 
doublets, but does not follow the above treatment, since supersymmetry ensures that 
there be no CP violation in the Higgs potential (A5 = Nonetheless, the MSSM has 
two possible sources of CP violation, the phases conventionally called 4>A an d 4>b, which 
are absent in the minimal standard model ||132| ; a combination of these will contribute to 
baryogenesis. The interactions in this model are given in terms of superfields by 

\U\uQH + DX D QH' + E\ e LH' + Id e~ i( ^ B HH'] _, 

_ _ ' (4.18) 

+ m 3 / 2 [\A\e^ A {UCuQH + D^ D QH' + E^eLH') + \fi B \HH'] A , 

along with mass terms for the gauginos. The CP violating phases 4>A and 4>b occur only 
in interactions involving the super-partners of the ordinary particles. 

Suggestions for how spontaneous baryogenesis in the MSSM at the weak phase tran- 
sition could proceed were first made in refs. [|51| and |52| ; the first quantitative analysis is 
given in ref. ]55|| . There it was found that baryogenesis was negligibly affected by the phase 
(pA, but proceeded due to the effects of the phase 4>b which appears in the mass matrix of 
the supersymmetric partners of the gauge bosons and the Higgs scalar s (the "mos"). 

The generalization of spontaneous baryogenesis from a single fermion to several 
fermions with a spacetime dependent mixing involves making a space-time dependent uni- 
tary change of basis on the fermions ipi — > Uij(x^)ipj in order to render the fermion mass 
terms everywhere real, positive and diagonal; the space-time dependence of U requires 
that we replace the kinetic energy terms in the Lagrangian by 

%.E. - £ K.E. + ^{rtidpUW . (4.19) 

It is simplest to choose U to be an anomaly free rotation, so that there are no WW terms 
generated, which forces (UHd^U) to couple to ordinary quarks as well as inos. In ref. ]55 
it is shown that the generalization of eq. ( [4.9p is 

3 [Mo u idU[(j>) 



T 



J CO #Trt/t(^)^^e k r a (^), (4.20) 



where (4>) co is the value of the Higgs vev for which F a starts falling exponentially (see eq. 
I)), and t^ k is the SU(2) Casimir matrix. Eq. (|4.20| ) leads to an acceptable value 



for n B /s for certain regions of parameter space in the MSSM. However, between the 
restrictions that the phase transition be sufficiently first order and that n B /s ~ 10 -10 , 
a large region of parameter space is excluded, and the theory is experimentally testable as 
a baryogenesis candidate, a luxury accorded to few models. 
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4-4- The nonadiabatic "thin wall" regime: the charge transport mechanism 

Thin bubble walls occur naturally in models with scalar fields with trilinear couplings 
(e.g., a familiar extension of the MSSM), as well as in the two Higgs model or singlet 
Majoron with some tuning of parameters. For bubble walls thinner than the strong in- 
teraction mean free path Q4.3j ), S w < It — 4/T, incoming fermions in the bubble wall 
frame interact with the wall like quantum mechanical particles scattering from a potential 
barrier. Due to CP violating interactions with the scalar field, the reflection coefficients 
for the fermions can be such that a net charge is reflected from the domain wall, which is 
transported into the region preceding the bubble. This charge asymmetry (not necessarily 
in baryon number) can bias anomalous 577(2) fluctuations in a large region in front of the 
bubble wall and produce a baryon asymmetry. This "charge transport mechanism" was 
introduced in the context of the singlet Majoron model, where the fermions scattering of 
the wall were massive neutrinos [ 48[ ; here we will summarize ref . ]53| on the nonadiabatic 



limit of the two Higgs model discussed above. Potentially the most compelling application 
of the nonadiabatic regime would be the MSSM extended by a sing let field PT| , but that 
model has not been analyzed to date. 

The procedure that has been developed for analyzing the thin wall scenario is as 
follows: 

(i) Compute the wall velocity and profile, including the space dependent CP violating 
phase, from the finite temperature effective potential. 

(ii) Calculate the reflection coefficients for fermions striking the wall with momentum k 
in the domain wall frame of reference by integrating the Dirac (Majorana) equation 
of motion. Because of CP violation, the reflection probabilities for particles and 
antiparticles can be different. 

(iii) Convolute the reflection probabilities with the incoming (boosted) thermal flux of par- 
ticles to determine the flux of fermions reflected from the wall. Due to CP violation, 
this flux will carry nonzero quantum numbers X. 

(iv) The nonzero X flux from the wall will lead to a nonzero charge density n x in the 
region preceding the wall. In general, a nonzero charge X will cause the free energy 
to be minimized for nonzero B, and one finds 



__ dF 



= -Cfe , (4-21) 
B=0 T 2 
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where nx is a nonzero AT-charge density. Given eq. (4TI), the baryon density in a 



region which is a distance d(t) = (z — t>u,t) in front of the domain wall equals 

9r a /■*/«« «9F 9r a /-^ , 

n B {z,t) = — — / dt ~dB =C ^T dtn x (z-v w t) 



oo 

OO 



9r a 1 

C — - — / dz n„(z) 
T v. ' v 



w JO 



(4.22) 



Thus one needs to calculate the total amount of charge in the region in front of the wall 
in order to compute the net baryon density produced. (As in the adiabatic scenario, 
baryons that are produced pass into the broken phase and become stable; a major 
difference though is that the baryon asymmetry is produced in front of the bubble in 
the symmetric phase, rather than within the bubble wall.) 



The applications of this procedure to the two Higgs model is detailed in ref. [p3| . Rather 
than compute the bubble wall parameters from the effective potential, we use a represen- 
tative wall profile 

1 + taxih(z/5 v 



e -iA0[tanh(z/5 w )-l]/2 23) 



where the top mass mt(T) = Afi>i and the wall width S w are taken as free parameters to 
vary, and A9 = —n is assumed; the resultant baryon asymmetry should scale linearly for 
small A9. The only fermion that has appreciable reflection probability from the domain 
wall is the top quark, and CP violation causes an asymmetry in the reflection rate between 
right-handed and left-handed top quarks, leading to a flux carrying net axial charge shown 
in fig. [5]. The asymmetry is maximized for particles whose longitudinal momentum in the 
wall frame is of order the top mass in the broken phase at that temperature (all particles get 
reflected or transmitted, if their longitudinal momenta are lower or much higher than mt, 
respectively) . From these k dependent reflection probabilities the momentum distribution 
of the reflected flux is computed. The transport of the axial top quark number in front 
of the wall is calculated by Monte Carlo techniques, allowing each reflected top quark to 
undergo Coulomb scattering via gluon exchange (with an infrared cutoff in the form of 
a plasmon mass) until the wall catches up. The proportionality factor C in eq. ( |4. 21|) is 
found to equal 4/(1 + 2n) in the standard model augmented to have n Higgs doublets, and 
from eq. (|4.22| ) we compute the baryon asymmetry. 

Axial top number is not a good quantum number, as it is violated by Higgs scattering, 



and so in ref. [|53| fermionic hypercharge was identified as the transported charge. In fig. 
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we display the results of the Monte Carlo simulation for several parameters, showing a 
substantial transport of charge in front of the domain wall. The penetration distance of the 
charge pushed in front of the wall depends on the velocity of the wall, and for v w = 1/ \f?>— 
the sound velocity — we find a distance of ~ 40/T; for slower walls, the distance is greater. 
There are several factors contributing to this large distance: (i) The incident particles 
selected by the reflection coefficients have longitudinal momentum pl ~ mt = 2T, and so 
have higher than average energy to start with; (ii) Upon reflection, the particle energies 
are enhanced by 7(1 + 0) — 2, for v w = l/y/3; (iii) Particles typically experience numerous 
collisions, but primarily in the forward direction. Our numerical results are consistent with 



Bjorken's energy loss formula ||133|| for light quarks in a quark-gluon plasma 



^ = ~a 2 a {l + N f /6)T 2 In ^/AET/M 2 , (4.24) 

where Nf = 6 is the number of quark flavors in the plasma and M is an infrared cuttoff (~ 
the plasmon mass): this formula yields a stopping distance of I ~ (9/87ro^)T -1 = 35T -1 
for a reflected top quark with initial energy E ~ 6T, twice the average thermal energy, 
(iv) Finally, for very slow domain walls, diffusion of charge in front of the wall plays a large 
role and increases the penetration distance substantially. 

Khlebikhov has shown by means of the Vlasov equations that the above scenario can- 
not proceed exactly as described above, since hypercharge is gauged and any hypercharge 
density produced will be screened by the plasma [|134|. In ref. [E3 this error was corrected. 



The correct charge to consider as the "X" charge in eq. (|4.22| ) is not fermionic hypercharge 
Yf, but rather the linear combination B' = B — xYf where B is baryon number and x is 
chosen so that B 1 is not screened by gauge forces; we find x = 2/(10 + n), where n is the 
number of scalar doublets. The reflected flux of top quarks carries B' charge, and upon 
redoing the analysis of ref. [53], the final results for the produced baryon asymmetry are 
found to be unchanged. These are shown in fig. [7[ The contour plot shows that a very large 
baryon asymmetry can result for maximal CP violation, v w = l/y/3 and 5 W x mt = 0{1), 
where mt is the value of the top mass in the broken phase at the phase transition. For a 



wall velocity v w = 0.1 one finds an order of magnitude enhancement in n B /s [[53]. 

We conclude that for thin bubble walls the charge transport mechanism works very well 
for producing the baryon asymmetry, and in fact suffers from fewer theoretical uncertainties 
than the adiabatic scenarios where baryogenesis actually takes place within the bubble wall. 
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4-5. Conclusions about mechanisms 

As is evident from the above discussion, the mechanism for weak scale baryogenesis 
depends on the width and velocity of the bubble wall. Allowing for certain processes to 
remain in thermal and chemical equilibrium during the phase transition, we have been able 
to show how weak scale baryogenesis could proceed in both the adiabatic and nonadiabatic 
regimes, producing a phenomenologically acceptable value for n B /s. It is tempting to "in- 
terpolate" and assume that baryogenesis would also work for the interesting case of phase 
transitions that lie between these two regimes, where scattering lengths are comparable 
to interaction lengths with the domain wall; however, the tools of equilibrium thermody- 
namics do not suffice for such a regime which would have to be analyzed by means of the 
Boltzman equation. Such a treatment could afford detailed knowledge about the crucial 
departure from equilibrium in all regimes. 

In addition to the models discussed above there are a number of others in the literature, 
which we include for completeness. An early model of McLerran [ 129 | attempted to use 
cosmological axions to provide the required CP violation in standard model baryogenesis, 
but the resultant baryon asymmetry proved to be too small. Other weak scale baryogenesis 
investigations include the effect of anomalous QCD events at the weak scale [|135| ; left- 
right models [ |136| , |137f |; transient scalar-top condensation in supersymmetric models [ 138 



lepton number production from domain walls colliding at the end of the transition [139 



and baryogenesis from electroweak strings without a first order phase transition [140 



5. CONSTRAINTS AND EXPERIMENTAL SIGNATURES 

One of the most appealing features of weak scale baryogenesis is the implication that 
new physics lies in wait for us at the weak scale. Since weak scale baryogenesis is feasible 
in most extensions of the standard model experimental guidance is needed. 

Possible and expected experimental consequences of weak scale baryogenesis are: 

1. New particles. The requirement that the weak phase transition be sufficiently first order, 
combined with Higgs mass limits from LEP, imply that the weak symmetry breaking sector 
be extended beyond a single Higgs doublet. Furthermore, the requirement that there be 
new CP violation also requires that new particles be added to the standard model. Both 
constraints suggest that there ought to be more than one scalar particle for the SSC to 
find. 
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2. New CP violation. The new sources of CP violation may be large enough to see in 
future experiments measuring the electric dipole moment of the electron and neutron [ |141| | , 
atomic physics ||142|| , or top quark production ||143|| . 

3. Constraints on supersymmetry. Spontaneous baryogenesis at the weak scale in the 
MSSM puts significant constraints on the experimentally allowed region of parameter 
space. For example the minimal supersymmetric standard model can avoid baryon number 
washout provided that the lightest Higgs scalar is lighter than about 50 GeV, tan/3, the 
ratio of the two Higgs vevs, is smaller than about 1.7, the top quark is heavier than about 
150 GeV, and the squarks are lighter than about 150 GeV ||106|| . (All of these bounds are 
only approximate because they are computed in one loop perturbation theory, and can be 
weakened by allowing the top squark to be significantly lighter than the other squarks.) 
Creating sufficient baryon asymmetry requires the electric dipole moment of the neutron 
to be greater than 10 -27 e-cm ]55 . 



4- Flavor changing neutral currents. Typical extensions of the standard model that are 
appropriate for baryogenesis (e.g. the two Higgs model) could have significant flavor chang- 
ing neutral currents large enough to see in the heavy flavor sector. Interesting places to 
look are B° -> jtx+jt - ||144|1 , t -> 07, and t -> cZ [p5j . 

5. Dark matter searches. Anomalous weak processes could explain the relative abundance 
of dark matter to baryon number in the universe today by coproduction of massive stable 
particles along with baryons through the anomaly ||146 [. These particles would have to 



carry weak interactions, and are tightly constrained by present dark matter searches [147 



6. FUTURE CALCULATIONS 

We have seen that the baryon number of the universe may have been produced at 
the weak phase transition via experimentally accessible weak scale physics. A precise 
calculation of the baryon abundance produced in a given model would allow us to pin 
down the parameters in that model, much as the abundances of the light elements and of 
dark matter are used to constrain models. So far, the baryogenesis calculations have been 
qualitative, and probably only accurate to within a couple of orders of magnitude. A more 
quantitative calculation will require further analytic and numerical work on the following 
questions: 
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1. A lattice study of the order and parameters (critical temperature, value of the order 
parameter at the critical temperature) of the electroweak phase transition. It is is of 
particular interest to see how the parameters depend on the Higgs mass. 

2. A computation of the baryon violation rate in both phases and within the bubble wall 
itself. 

3. A quantitative calculation of the nonequilibrium transport of quantum numbers is 
needed both for all but the most adiabatic cases, for instance by numerically solving 
the Boltzmann equations. 
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Figure Captions 



Fig. 1. The potential energy of the 577(2) gauge field along the "sphaleron" trajectory. 

The minima correspond to configurations with zero gauge field strength, but 
different baryon number. The lower curve ignores the energy of the baryons, 
while the upper curve does not. The zero of energy for the two curves has been 
shifted for clarity. 

Fig. 2. Ncs as a function of time in the Abelian Higgs model. (From ||62|| ). 

Fig. 3. The effective potential for a theory which has a first order phase transition. The 
dotted line is the potential at the critical temperature T c , which has two degen- 
erate minima. The solid line is the potential at the low temperature To, where 
only the broken symmetry phase is stable. The dashed line shows the potential 
at an intermediate temperature where the broken phase has lowest free energy 
but the symmetric phase is metastable. The dot-dashed line is the potential at 
high temperatures, where only the symmetric phase is stable. 

Fig. 4. Typical "daisy" contribution to the effective potential 

Fig. 5. The difference in reflection probabilities for left and right handed top quarks, as 
a function of the longitudinal momentum in the wall frame, for three different 
wall widths. The wall shape is given in ( j4.23| ), with A9 = — n, corresponding to 
maximal CP violation. 

Fig. 6. A plot of the ratio of hypercharge to entropy densities preceding the bubble wall, 
as seen in the thermal frame of reference. The horizontal axis is the distance 
ahead of the wall, in units of the inverse temperature. These computations were 
performed assuming rrit = 2T and 5 W = 1/T, for the two sample wall velocities 
v w = 1/a/3 (7 = 1.225) and v w = .98 (7 = 5.0). It is the area under these curves 
that enters eq. ( f4.22|) for n B . 

Fig. 7. Contour plot of n B /s as a function of the top mass and 5 W , in units of the critical 
temperature, assuming maximal CP violation (A9 = — it in eq. ( 4.23|) ). We 



have assumed the number of light scalars to be n = 2, the wall velocity to be 
v w = l/y/3 and have taken k = 1 in eq. (|2.9| ). 
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